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Abstract 

We present calculations for the temperature-dependent electronic structure of bulk 
ferromagnetic EuO based on a parametrization of the d-f model Hamiltonian us- 
ing results of first-principles TB-LMTO band structure calculations. The presented 
method avoids the problem of double-counting of relevant interactions and takes 
into account the symmetry of the atomic orbitals. It enables us to determine the 
temperature-dependent band structure of EuO over the entire temperature range. 
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The model of choice for describing the electronic and magnetic properties of 
magnetic rare-earth systems is the s-/(or s-d) model [1]. The model describes 
the exchange coupling of itinerant band electrons to localized magnetic mo- 
ments. It is applicable to the magnetic semiconductors as the europium chalco- 
genides and to the metallic elemental rare earths Gd, Tb, and Dy. Many 
characteristics of these materials can be explained by a correlation between 
the localized magnetic -^/states and extended conduction band states, mainly 
5d-6s. Within the s-f model, the correlation is represented by an intraatomic 
exchange interaction. 

There have been earlier works on the temperature-dependent band structure 
of bulk EuO [2,3]. The presented approach employs a decomposition of the 
relevant Eu-5d bands of EuO into five consecutive bands. For each of these 
bands, a single-band s-/ model calculation has been performed. By employing 
this splitting procedure, the method employed by the authors [2,3] disregards 
the full multi-band structure of the conduction bands in EuO and leads to an 
overestimation of relevant correlation effects. 

In the following, we will introduce a multi-band d-/ model and use it to calcu- 
late the temperature-dependent bandstructure of EuO. The employed method 
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combines a reliable many-body treatment of the d-f exchange model with a 
self-consistent band-structure calculation. In particular, the full multi-band 
aspect of the problem is treated in a reliable and symmetry-conserving mode. 



1 Theory 

The conventional or single-band s-/ model describes the intratomic exchange 
interaction between an s-like conduction band and system of localized mag- 
netic moments, which are due to the partially filled .^/-shells. For a real ma- 
terial, like EuO, the conduction band usually consists of multiple subbands, 
which can be distinguished by their magnetic quantum number. In the follow- 
ing we will introduce the model for such a local-moment system with multiple 
conduction bands. The respective model will be referred to as the (multi-band) 
d-f model as opposed to the (single-band) s-f model. 

The d-/ model Hamiltonian consists of three parts: 

H^Hd + Hf + Hdf. (1) 

For a real system like EuO, the kinetic energy of the conduction electrons 
reads: 

H,= Y: Tr'4maC3m'a. (2) 
ijmm' 

Here, cl^„ and Cima are, respectively, the creation and annihilation operators 
of an electron from the m-th subband at the lattice site Rj. T^'^' are the 
hopping integrals. 

Each lattice site Rj is occupied by a localized magnetic moment, represented 
by a spin operator Sj. These localized moments are exchange coupled expressed 
by the Heisenberg Hamiltonian: 

= Y JijSi ■ Sj, (3) 

ij 

where Jij are the Heisenberg exchange integrals. 

In addition to the contribution of the itinerant-electron system and the contri- 
bution of the localized /-moments we have a third contribution to account for 
an intraatomic interaction between the conduction electrons and the localized 
/■spins. The form of that contribution can be derived from the general form 
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of the CoTilomb interaction between electrons at a single lattice sites. If the 
electron scattering processes caused by the Coulomb interaction are restricted 
to two involved subbands L and L', one gets for the Coulomb interaction: 

Hl^\Y^ [ULL'ctaCL'a'CL'cr'CLc7 + JLL'ctaCt'a'CLa'CL'a 
LL'aa' 

+ JlL'CtaCL'CL'a'CL'a]- (4) 

We can split the Coulomb interactions into three different parts, depending on 
whether both L and L' are conduction bands, Hf^, both belong to localized 
bands, Hj'^, or there is an interaction between localized and conduction bands 
involved, Hf^, 

Hi = Hf^ + H^^ + Hf^. (5) 



The first of the three parts disappears, since for the case of a semiconductor 
the conduction bands are unoccupied. The second term is already contained in 
the description of the locahzed moments via the Heisenberg model in Eq. (3). 
The third term is the only one which contains an interaction between the 
localized bands and the conduction bands. Assigning the indices m and / to 
itinerant and localized bands, respectively, it can be written in the form: 



Hi^ = X] [ UmfC^aCfa'Cfa'Cma + JmfC^aCfa'Cma'Cfa 
mfaa' 

The last two terms of the above equation vanish for the special case of EuO, 
due to the two fact that in Eu the localized 4f shell has its maximum spin of 
S = 7/2 resulting in no double occupancy of the different subbands of the 4f 
shell. Introducing the Pauli spin operators ctl one arrives at 

Hf^ ^ —2^Jjnfcrjn - af. (7) 

mf 



By defining the spin operator of the local moment by S = fiY^fCTf and by 
assuming the inter-band exchange to be independent on the band indices m 
and /, Jmf = J/2, the multi-band (/-/interaction which is the sum of Eq. (7) 
over all the different lattice sites reads: 

Hdf = -T- X] '^irn • Sj. (8) 
im 

where J is the intraatomic c?-/ exchange interaction. Using the second-quantized 
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form of (7j and the abbreviations 



Sj — Sj iZ(j-Sj'^ 



(9) 



the d-/ Hamiltonian can be written as 




(10) 



The most decisive part of the d-/ Hamiltonian (10) is the second term, which 
describes spin exchange processes between the conduction electrons (2) and 
the locahzed moments (3). 

The many-body problem that arrises with the Hamiltonian (1) is far from 
being trivial and a full solution is lacking. In previous papers we have presented 
an approximate treatment for single-band s-/ model in a local- moment system 
with film geometry for the special case of an empty conduction band, n = 
[4]. For this special case, the approach for the single-band s-f model in a 
film geometry with inequivalent layers a can be identically transferred to the 
situation of the multi-band d-f model in bulk materials with different subbands 
m, via the index transition a ^ m. Using this transition, in the following we 
want to present the main results of the theoretical approach for the multi-band 
d-/ model. 

Due to the empty conduction bands we are considering throughout the whole 
paper, the Hamiltonian (1) can be split into an electronic part, Hd + Hdf, and 
a magnetic part, Hf, which can be solved separately [4]. For the magnetic sub- 
system Hf we employ an approach based on the random phase approximation 
(RPA) which has been described in detail in [5]. As the result one gets the 
temperature-dependent magnetizations of the local-moment system. 

Concerning the electronic subsystem, + Haj, all physically relevant infor- 
mation can be derived from the retarded single-electron Green function 



The Fourier transform of the sing le-electron Green function, (E) = iV"! Eii e^'^(^'-^^)G^™'(E), 
is related to the spectral density via 




E 



oo 




(11) 







(12) 
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from which the partial density of states of the m-th subband can be obtained: 

Pa{E) = ^ESnE). (13) 

For the solution of the many-body problem posed by the Hamiltonian Hd+H^f 
we write down the equation of motion of the single-electron Green function 

(11) 

E G^f = + E ^r"G^^' + (([q^., HafU cU'Je. (14) 

km" 



The formal solution of Eq. (14) can be found by introducing the self-energy 

"[Cima. H,f]_ ; c^,^))^ = E Ml^f{E)Glff{E), (15) 

km" 



which contains all information about the correlations between the conduction 
band and localized moments. After combining Eqs. (14) and (15) and perform- 
ing a Fourier transform we see that the formal solution of Eq. (14) is given 

by 

Gka(£;) ^h{EI-T^- M^,{E)y^ , (16) 



where I represents the (M x M) identity matrix and where M is the number 
of subbands of the condcution-band system and where the matrices Gk„{E), 
Tk, and Mko-(-E') have as elements the subband-dependent functions {E), 
T^"^', and M{^J^'{E), respectively 

The evaluation of the higher Green functions resulting from (14) employs a 
moment- conserving decoupling approximation (MGDA). Following the steps 
presented in [4] one arrives at an implicite set of equations for the local self- 
energy, M^J-' (E) = 5mm'M^{E): 

Mr(^) = -^//-(£;), /.n^) = ^|||, (17) 



where the numerator and the denominator, respectively, are given by 



Z^^z,h^{S^) + |{(4^) - h^S{S-rl))G\ 



mm 
Oct 



2 

((aw + A(^) + z^^iZ)Hs^) + ^4'^)Gr.} 
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+ -^[z„h'S{S + 1)(A« + A(2) + z./x-J 

+Kf (/.^ - ^i'!l,)}G^G^a, (18) 

+zAf^ i^a + fi)}GrGr, (19) 

where 



K - ^z^^ , (20) 

contain /-spin correlation functions, which have to be calculated within the 
Heisenberg model. 



For the limiting case of ferromagnetic saturation of the localized ^spin system 
{T — 0), there exists an exact solution for the d-f (s-f) model with empty 
condcution band (n = 0) [6]. In particular, the spectrum of the spin-| electron 
is only rigidly shifted compared to the free solution obtained for the case of 
vanishing s-/ interaction, J = 0. This means, that the spin-| spectra obtained 
within a ferromagnetic band-structure calculation provide the single-particle 
input (cf. Eq. (2)) for the d-/ model calculation. Since the solution of the d-f 
for T = and n = is exact, the problem of double counting of relevant 
interactions is elegantly avoided in the presented approach. 

To obtain the temperature-dependent band structure of EuO, the T — 
band structures are needed as single-particle input for the model calculations. 
These have been calculated for the Eu- 5d bands, which in EuO dominate the 
conduction band, using a standard TB-LMTO program [7]. The results of 
these calculations will be presented in the next section. 



2 Band-structure calculations 



The europium chalcogenides crystallize in the fee rock-salt structure. For EuO, 
the lattice constant is a = 5.142 A = 9.717 au [8,9]. The electronic confi dura- 
tions of europium and oxygen in the state of the free atom and in the ionic 
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picture for EuO, respectively, are 

free: [Eu] = [Xe]{4fy{Qsf, [O] = [Be]{2sf{2pf; 
EuO: [Eu] = [Xe](4/)^ [O] = [Ne]. 

In accordance with Hund's coupling, the Eu-^/ electrons couple to the maxi- 
mum magnetic moment of 5 = 7/2. The difficulty in dealing with the 4f levels 
within an LDA calculation lies in their strongly localized character, which is 
due to a strong Coulomb interaction between the 4f electrons. As a result 
of the inability of the LDA to take into account the respective interactions 
correctly, a normal LDA calculation for EuO produces a metal with the 4f 
levels lying well within the conduction band. 

To overcome this situation, we chose to treat the 4f moments as localized core 
electrons^. Fig. 1 shows the respective spin-dependent band structure of EuO. 
In Fig. 2, for the same calculation, the density of states is displayed. Here, the 
main contributions to the density of states originate from the 0-2p levels, 
which constitute the valence band, and from the unoccupied Eu- 5d levels. As 
a result of treating the 4f levels as core electrons, we have the correct ground 
state of an insulator. Clearly, the conduction-band region is dominated by the 
Eu-5d electrons. Thus, for the model calculations which will be performed 
in Sec. 3 it is a reasonable approximation to restrict the single-particle input 
obtained within the band-structure calculations to the Eu-5c? bands. 



3 Temperature-dependent band structures 

In this section the LSDA band structure of bulk EuO shall be combined with 
an d-f model calculation to obtain the temperature-dependent band struc- 
ture. Due to the fact that within our model calculation for T = the spin-| 
spectrum is rigidly shifted towards lower energies by the amount of ^JS, the 
obtained k-dependent hopping matrices for the spin-| electron will serve us as 
the input for the kinetic Hamiltonian (2). Since the solution for T = is exact, 
the problem of double counting of relevant interactions, which usually occurs 
when combining first-principles and model calculations, is elegantly avoided. 

The temperature comes into play via the temperature-dependent /spin corre- 
lation functions of Eqs. (20), which can be calculated within the Heisenberg 

^ For EuO, it has been shown [10] that the band structures of the relevant conduc- 
tion bands obtained within such a calculation differ only slightly from those obtained 
within an LDA+U caculation. For the 4f system Gd the different approaches for 
dealing with the localized 4f moments have been reviewed by Eriksson et a/. [11]. 



7 




Fig. 1. Spin-dependent (spin-|: , spin-|: - - -) band structure of bulk EuO 

calculated within an LSDA calculation with the 4f levels treated as core electrons. 
The horizontal dashed line represents the Fermi energy. 

model for the local-moment system [10,4,5]. For the following investigations, 
however, the concrete shape of the magnetization curve is of lesser interest 
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Fig. 2. Density of states, Pa{E), of bulk EuO calculated within an LSDA calculation 
with the 4f levels treated as core electrons. The upper graph shows the total density 

of states, whereas the lower two display the partial density of states of the 0-2p 
and the Eu-5d bands. The dotted line represents the Fermi energy. The densities of 
states above the Fermi energy are multiplied by a factor of three and correspond to 
the left y-axis. 

and the reduced magnetization {Sz)/S has been chosen as the temperature 
parameter. 

To calculate the temperature-dependent band structure of EuO, the d-/ ex- 
change interaction is needed. It has been shown by several authors[12-14], that 
the standard LSDA band-structure calculations are quite compatible with 
the simple Stoner or mean-field picture, in which the exchange splitting is 
only slightly energy dependent. Contrary to this statement. Fig. 2 implies an 
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energy-dependent splitting of the spectra. However, the inclusion of a such an 
energy-dependent splitting of the spectra is not possible within the theory pre- 
sented in this paper and we have to take an averaged exchange interaction for 
the calculation of the temperature-dependent band structures. Averaging the 
splitting of the spectra for the two spin directions over the energy, one arrives 
at an average splitting of the spectra for the two spin directions of 0.875 eV. 
In the mean-field approximation of the d-f model for T = the spin-splitting 
is equivalent to JS. With S = 7/2 one gets for the d-f exchange splitting 
J = 0.25 eV. This value will be used for the further d-f model calculations. 

Fig. 3 shows the temperature-dependent spectral density of bulk EuO obtained 
within an d-f model calculation with J = 0.25 eV. For T = {{S,)/S = 1), 
the spin-| spectral density agrees, except for a constant energy shift, with 
that obtained by the LSDA calculation. However, still for T = but for the 
spin- J, spectrum a broadening of the dispersion curves can be observed, most 
notably around the F-point, which indicates a finite lifetime of the respective 
quasiparticles due to correlation effects. In this respect already the T — 
spin- 1 solution goes beyond LSDA by taking into account correlation more 
realistically. 

For intermediate temperatures {{Sz)/S = 0.5) a broadening of the dispersion 
curves sets in also for the spin-| spectra since the spin-| electron can now ex- 
change its spin with the deviated local- moment system. For the spin- J, spectral 
density, the correlation effects, already present at T = 0, increase, as can be 
seen by the further broadening of the curves. At the same time, the spin-| and 
the spin- 1 spectra are shifted towards higher and lower energies, respectively, 
reducing the effective splitting between the two spectra. This effect contin- 
ues with increasing temperature, until finally, for T — Tq {{Sz)/S — 0), the 
spectra for both spin directions are equal. Clearly, the temperature-dependent 
effects displayed in Fig. 3 do not comply with the simple Stoner-picture which 
would imply a constant energy shift of the spectra. The reason for this more 
complicated behavior as a function of temperature is again that in the theory, 
the correlation is treated in a way which goes beyond mean-field. 

The mentioned shift of the spin-| spectrum towards lower energies when going 
from T = Tq down to T = represents the red shift of the optical absorption 
edge in EuO [15, 16]. Fig. 4 displays the densities of states obtained for different 
magnetizations of the 4f moments. It can be seen that the red shift of the 
optical absorption edge is due to the 4f-5dt2g transition. From Fig. 4 one 
obtains a red shift of rs = 0.35 eV. This value agrees reasonably with the 
experimental value for the red shift of 0.27 eV [9]. Here, the agreement can be 
further improved when choosing the splitting of the lower edges of the LSDA 
Eu- 5(i bands in Fig. 2 to calculate the d-f exchange interaction for the model 
calculation. 
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Fig. 3. Spin-dependent spectral densities of the Eu-Jd bands of bulk EuO for 
J = 0.25 eV and for different magnetizations {Sz)/S. 



There have been previous calculations by Nolting et al. concerning the temperature- 
dependent band structure of bulk EuO [2,3] . In these works, the band structure 
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Fig. 4. Temperature-dependent densities of states of the Eu-5d bands of bulk EuO. 
For r = {{Sz)/S = 1) the spectra for the two spin directions are furthest away 
from each other and approach each other when increasing the temperature. For 
T = Tc {{Sz)/S = 0) the densities of states for both spin directions are the same 
(fat lines). 

of the Eu- 5d bands has been split into five s-fike bands with the lowest eigen- 
states belonging to the lowest band etc. The splitting into 5-bands produces 
five consecutive bands of an average bandwidth of about = 3eV[2,3]. 

Contrary to these works, for the calculation of Figs. 3 and 4 the full band 
structure of the Eu-5d bands has been taken into account, thereby respecting 
the symmetry of the different Eu-Sd orbitals. As a result, the Eu-5d bands of 
bulk EuO split into the t2g and the Cg subbands with a bandwidth of about 
10 eV and 6eV, respectively. For the model calculations, the decisive entity 
for the magnitude of the correlation effects is the d-f exchange interaction 
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over bandwidth, J /W . In this respect, in the prcvioTis works [2,3] the calcu- 
lated correlation effects should be slightly overestimated, considering the more 
"natural" band decomposition employed in this thesis. 



4 Conclusion 

We presented a method for calculating the temperature-dependent band-structure 
of a local-moment system based on a parametrization of the (multi-band) d-f 
model Hamiltonian using results of first-principles TB-LMTO band structure 
calculations. The method has been applied to the ferromagnetic semiconductor 
EuO. 

The presented method enables us to calculate the temperature-dependent 
bandstructure of EuO over the entire temperature range. A main characteris- 
tics of the approach is the avoidance of double-counting of relevant interactions 
for the conbination of first-principles and model-calculation. 

Concerning the earlier works on the temperature-dependent band structure of 
bulk EuO [2,3], the method presented in this paper has two main advantages: 

(1) The full multi-band character of the conduction bands, including in par- 
ticular the symmetry of the relevant orbitals, is taken into account. 

(2) Due to the broader bandwidths of the involved subbands the correlation 
effects should be treated in a more realistic fashion. 

The extension of the method to systems with finite band occupation, n 7^ 0, 
and to systems with reduced dimensionality promises to give an insight into 
the fascinating physics at the surface of rare-earth systems [17]. 
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